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GENERIC PROPERTIES FOR RANDOM REPEATED QUANTUM 

ITERATIONS 


ARTUR O. LOPES AND MARCOS SEBASTIANI 


ABSTRACT. We denote by M" the set of n by n complex matrices. Given a fixed density 
matrix j3 : C" —» C' 1 and a fixed unitary operator U : C" <g> C" — » C" 0 C", the transformation 
& :M" —> M" 

Q^4>{Q) = Tr 2 {U(Q®p)U*) 

describes the interaction of Q with the external source /3. The result of this is &(Q). If Q is 
a density operator then is also a density operator. The main interest is to know what 
happen when we repeat several times the action of in an initial fixed density operator Qq. 
This procedure is known as random repeated quantum iterations and is of course related to 
the existence of one or more fixed points for <E>. 

In 0, among other things, the authors show that for a fixed there exists a set of full 
probability for the Haar measure such that the unitary operator U satisfies the property that 
for the associated d> there is a unique fixed point Q$>. Moreover, there exists convergence 
of the iterates <3> n (<2o) —> Q<t>, when n —> °o, for any given Qq 

We show here that there is an open and dense set of unitary operators U : C n <S> C n —> 
C n ® C" such that the associated has a unique fixed point. 

We will also consider a detailed analysis of the case when n — 2. We will be able to 
show explicit results. We consider the C° topology on the coefficients of U. In this case we 
will exhibit the explicit expression on the coefficients of U which assures the existence of 
a unique fixed point for Moreover, we present the explicit expression of the fixed point 
(2<t> 


1. Introduction 

We denote by M n the set of n by n complex matrices. Given a fixed density matrix 
/3 : C” —>• C” and a fixed unitary operator U : C n 0 CM C w 0 C", the transformation 
: M n -» M n 

Q^<P(Q) = Tr 2 (U(Q®ll)U*) 

describes the interaction of Q with the external source j3. 

We assume that all eigenvalues of p are strictly positive. 

In 01 the model is precisely explained: Q is in the small system and /3 describes the 
environment. Then <f>((7) gives the output of the action of p m Q given the action of the 
unitary operator U. 

Other related papers are 12) and 0. 

The main question is about the convergence of the iterates 0" ( Qo ). when n —> for 

any given Qq. It is natural to expect that any limit (if exists) is a fixed point for 0. 

Our purpose is to show the following theorem: 

Theorem 1. Given a fixed density matrix p : C" —> C n ,for an open and dense set of unitary 
operators U : C" ® C" —> C" (g) C" the transformation 0 : M n —y M n 

Q^<$>{Q)=Tr 2 (U(Q®p)U*) 


l 
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has a unique fixed point Q<t>. In the case n = 2we present explicitly the analytic character¬ 
ization of such family ofU and also the explicit formula for Q&. 


This result implies one of the main results in 0 that we mentioned before. 


2. The general dimensional case 

Suppose V is a complex Hilbert space of dimension n > 2 and (V) denotes the space 

of linear transformations of V in itself. 

Then, Tr 2 : J£?(V®V) Jz?(V), given by Tr 2 (A®B) = Tr(B)A. 

There is a canonical way to extend the inner product on V to V 0V. 

We fix a density matrix /3 £ 2zf(V). For each unitary operator U £ 2z?(V 0 V) we denote 
by T>[/ : Jf(V) —> JSf (V) the linear transformation 

<Pu(A) = Tr 2 (U(®P)U*). 

We denote by T C Jz? (V) the set of density operators. It will be shown that <tv preserves 
r. As r is a convex compact space it has a fixed point. 

The set of unitary operators is denoted by f/. 

If A is such that <P/y (A) =A, then it follows that the range of O/j — / is smaller or equal 
to n 2 — 1. 

We will show that there exist a proper real analytic subset Xcf such that if U is not 
in X, then range <J>[/ — I = n 2 — 1. In this case the fixed point is unique. More precisely 

X = {U £ f/ : range (Oy —/) < n 2 - 1}. 

This X C is an analytic set because is described by equations given by the deter¬ 
minant of minors equal to zero. It is known that the complement of an analytic set, also 
known as a Zariski open set, is empty or is open and dense on the analytic manifold (see 
ID). Therefore, in order to prove our main result we have to present an explicit U such that 
range of (<tv — I) is n 2 — 1. 

This will be the purpose of our reasoning described below. 

The bilinear transformation (A,B)—>T r(B)A from ££ (V) x 2z? (V) to 2z? (V) induces the 
linear transformation 


Denote by e\,e 2 , ...,e„ an orthonormal basis for V. We also denote Ljj € the 

transformation such that Li fief) = e,- and Ljjfifi = 0 if kfi j. 

The Lij provides a basis for.Z'(y). 

If A £ jSf(y) we can write A = a//T,y and we call [afi \ <ij< n the matrix of A. 

Note that e,- <g) 1 < i. j < n is an orthonormal basis of V (dV . Moreover, 

Lik®Lj,{e k ®ei) =ei®ej, 

and 

Lik®Lji(e p ®e q ) = 0 if (p,q) ( k,l ). 

It is also true that: 

a) Ljj Lpq — 0 if / f p, 

b) Li j Lpj — Liq , 

c) Tr (Lij) = 0 if if j and Tr (La) = 1 . 

One can see that L,*. ®Lji, 1 < i,k,j,l < n is a basis for fif(V <8> V). 
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Given T £ JS?(V Cg> V) denote T = Lhjft ; /^ft ®Lji. Then, 

Tr 2 (T) = = £ (L'.jftj)**- 

ft i 

In the appendix we give a direct proof that: if A £ T, then 4>y (A) £ r, for all U £ . 

Now we will express <tv in coordinates. We choose an orthonormal base e\.e 2 .... e n £ V 
which diagonalize /3. That is 

J3 = Y \ L m- \ >0, 1 < q < n, Y\ = 1. 

9 9 

Given r,s, 1 < r. s < n, we will calculate <t>u(L rs ). 

Suppose U = 'L u i.j,k-l^ik ®Lji, then U* = Y. Li ij,k,iL ik ®Lj t and 

(J-'rs 0/3)1/ = (Yj Lrs ® Lqq') U = U k ,l,sJ L rk ® 

9 j 

Now, we write U = Y* u a,p,Y,dLay ®Lpg. Then, we get 

U(L rs 0 /?) t/ u k,l,sjLak ® Lpj. 

Finally, 

^u(^rs) = Yi^J Ua Ji r J U k,l,s,j L a k = Y* ( Yi^j Ufx J> r J U k,l,s,j ) ^oft- 

aft j',/ 

As r is convex and compact and (f>u is continuous as we said before there exist a fixed 
point A £ r. In particular the range of (j)u is smaller or equal to n 2 — 1 . 

We will present an explicit U such that range of (<S>j/ — I) is n 2 — 1 . 

This will be described by a certain kind of circulant unitary operator 

Suppose mi , M 2 , • • •, u n 2 are complex number of modulus 1 . We define U in the following 
way 

U(e i ®ei) = mi (ei 0e2), U(e i ®e 2 ) = U 2 (<?i ® e 2 ), ...,U{e\ ® e„) = u„ (e 2 ®ei), 


U(e 2 ®ei) =u n+ i(e 2 ®e 2 ), U(e 2 ®e 2 ) = u n+ 2 (e 2 ®e 2 ),...,U{e 2 ®e n ) = u 2 n (e 2 ®e i), 


U(e n ®e\) = u n 2_ n+ 1 {e„®e 2 ), U(e n ®e 2 ) = u n 2_ n+ 2 (e n ®e 2 ,),...,U{e n ®e n ) =u „ 2 (e\®ei), 

We will show that for some convenient choice of u \, u 2 ,u n i we will get that the range 
of <!>[/ —/is n 2 — 1 . 

Suppose 


G — U’i.j.kJ Ljk V) Lji . 

in this case 

U(e k ®ei)=Y* u iJM e i ® e j- 

ij 

By definition of U we get 

a) if l < n, then m,- ,■ */ ^ 0 , if and only if, i = k, j = l + 1 ; 

b) if k < n, then Uij.k,n ^ 0 , if and only if, i = k+ 1 , j = 1 ; 

c) Ui,j,n,n 7 ^ 0 , if and only if, i = j = 1 . 
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For fixed r,s such that 1 <,r,s <n we get from a), b) and c): 

1 < r < n, 1 < s < n, implies 

n — 1 

&U(Lrs) ~ ( ^ I .r.j ^.v,/+1 ,s.j Ay ) Lrs T Mr l.l j'.n M$ j-1.1 ,s.n 'Kn ^(r+l)(s+l); 

7=1 

1 <s<n, implies 

n— 1 

^ u(Lns ) = ( ^2 ^n.j I .n.j ^s,j+l,sj 4/ ) Lns "1” lfi.1 ,m,m 1.1 ,5,n An ^1 ( 5 + 1)5 

7=1 

1 < r < n, implies 

n— 1 

3>£/(L„) = ( ^2 llr.y+1 .r.y U n J+l,nJ Ay )An T M r+l,l,r,n Ml,l,n,n A n ^(?'+l)l- 

7=1 

In particular for 1 < r < « we have <I>y (L, T ) = ( 1 — X n )L rr + A„L( r+ i j ( r+ i). 

In order to show that the range of < 3 >j/ — / is n 2 — 1 we will show that the (j>u (L rs ) — L rs 
are linearly independent for ( r,s ) ^ (n,n) 

Suppose that 

^ c„(fo/(L„)-L„)= 0 . 

(r,5)^(n,n) 

The coefficient of Ln is —A„cn, then ci i = 0 . 

The coefficient of L22 is A„cn — A„C22, then C22 = 0 . 

The coefficient of L nn is A„C(„_i)( n _i), then C( ;! _n(„-i) = 0 . 

Then, we get that 

Y,Crs(<l>u{Lrs)-Lrs) = 0 . ( 1 ) 

r^s 

We will divide the proof in several different cases, 
a) Case n = 2 . 

52 c rs (A (Cm) — Cm) = C12 ( 0 f/(Cl 2 ) — C12) +C21 (•)>[/(C21) — C21). 

r^s 

By definition of £/ we have that u i 2, 1,1 = u 1, W2, 1,1,2 = ^2 5 ^2.2,2.1 = ^3 ? ^1,1,2,2 = ^4 • 
Therefore, 


and 


0t/ (C12) C[2 = (m 1 M3 Ai — l)Li2 + M 2 M 4 A 2 C 2 i 


<t>u(L 2l) — C21 = (M3 Ml Ai — 1)^21 + M4 M2 A2 Ci2- 
From ([!} it follows that 


(m 1 M3 Ai — 1) C12 + M4M2 A2C21 = 0 


M 2 M 4 A 2 C 12 + (M 3 mi Ai — 1 ) C 21 = 0 . 

Taking U such that it \ = i, 112 = M 3 = M 4 = 1 it is easy to see that the determinant of the 
above system is not equal to zero. Then we get that £+2 = C 21 = 0. 

Then, we get a U with maximal range. 
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b) Case n > 2. 

We choose ui,U2,---,u n 2 according to LemmaQ]below. 

The equations we consider before can be written as 

1 1 ^ s < n, r ^ s, then, Of/ {Lrs') L rs = {u rs 1) L rs T (s+i)) 

1 < s <C n, then, Oj/ ( L ns ) L ns — {ci ns 1) L ns -F b ns L |^^, 

1 then, Of/ {L rn ) = {u rn 1) H- b rn L^ r ^.\'j \. 

For instance 


and 


H-l 

®rs ~ ^r,j+l,rj Ms,j+l,s,j X j : 
7=1 


brs — 

Note that u r j + \ r j u s .j+\ }S j has modulus one and also u r +i.\. r .n «s+i.i,s,«- 
Moreover, | b rs \ = X„ > 0 and | a rs | < Ai +... + X n \. Indeed, note first that the products 
u r j+i.r.j u s j+i,s.j are different by the choice of the Ui ; k i (see LemmaQ}. Furthermore, by 
Lemma[2]we get that |a rs | can not be equal to Ai +... + A„_i. 

Therefore, | a rs — 11 > 1 — |a r j| > 1 — X q = A„ = | bij | > 0, for all r, s, i, j and r^s, 

i ± j- 


Suppose 2 < k < n. 

Remember that the L, ; define a linear independent set. 

The coefficient of L\ k in i[T} is 

Cu(«U~ 1) + C n (k-l)bn{k-l) = 0- 

The coefficient of i) in 01 is 

c n(k-l) ( a n(k-l ) — 1) + c («-l) (t-2) ^(n-1) (i-1) = 0- 

The coefficient of L^ n _ k+1 ) j in <|T} is 

C{n—k+2) 1 {®(n—k+ 2) 1 1) T C(n—k— 1)h^(h— k+\)n 0. 

The coefficient of Lt n _ k+ 1 )„ in £[} is 

C(n—k+\)n ( a (n—k+l)n ( ) X~ C(n—k)(n— l)^(n— k)(n— 1) 0- 

The coefficient of £(„_/•)(„- 1 ) in <[T]) is 

C{n—k) (n— 1) (^(n— k) (n— 1) () A" G(n—1) (n—2) ^{rt—k— 1) (n—2) 0* 


The coefficient of A> 2 (*+i) i n 01 is 

c 2(/fc+l) ( a 2(k+l) — 1) + c\ k b\ k = 0. 

If c u ± 0, then, from above we get |c u | < \c n{k _i)\ < ... < \c 2 (k+i)\ < ku|- 
Then, we get a contradiction. It follows that cj k = 0. 

Therefore, 

Cn(k— 1) H(n—1) (k~2) ^(w— k+2) 1 ^fn—fc+l)w ^(n—k) (n— 1) ••• ^2{k+\) 

From this follows that c rs = 0 for all r. s, when r / .v. This shows that for such [/ we 
have maximal range equal to n 2 — 1. 

Now we will prove two Lemmas that we used before. 
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Lemma 1. Given m > 2, there exists complex numbers u\..... u m of modulus 1, such that, 
if l < i j < m, 1 < k l < m and uiU] = u^Tf, then i = k,j = l. 

Proof: The proof is by induction on m 
For m = 2, just take u\ U 2 not in M. 

Suppose the claim is true for m > 2 and u\ , the corresponding ones. 

Consider 

S = {uiTTj | 1 < i,j < m } 

and 

T = {upUq | 1 < p,q < m }. 

Then, take u m+ \ such that u m+ iTTjl is not in S for all 1 < p < m, and n^ +1 is not in T. 
Then, u\, ...,u m ,u m+ \ satisfies the claim. 

□ 

Lemma 2. Consider A],.... A,„, real positive numbers and z \,.complex numbers of 
modulus 1. 

Suppose | Y!j= \ ^jZj I = Y!j= i kj, then zi = Z 2 = = Z m - 

Proof: The proof is by induction on m. 

It is obviously true for m = 1 . 

Suppose the claim is true for m — 1 and we will show is true for m. 

Note that 

m m m— 1 m 

L X > = \Lla l<l E l+<IV 

j =i j =i j =i j =i 

From, this follows that 

m— 1 m— 1 

11 ^ 1=1 Xj. 

7=1 7=1 

Then, zi = z 2 = ... = z m -1 = z. 

Therefore, 

m m— 1 m —1 m 

= Zm Kn | < |Z ^ A/ | + |z m A m | = ^ Ay. 

7=1 7=1 7=1 7=1 

Given t'i,V 2 complex numbers such that |vi +V 2 I = |vi| + |V 2 I, then they have the same 
argument. 

Then, there exists an s > 0 such that z A j = .vz m A m . 

Now, taking modulus in both sides of the expression above, we get 

m ~1 m—1 

A j — | 7 , Ay | = .V Z/7/ fni = sA m . 

7=1 7=1 


From this follows that z m —Z 


□ 
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3. The two dimensional case - explicit results 


Our main interest in this section is the explicit expression of the U such that the fixed 
point is unique. We restrict ourselves to the two dimensional case. 

We will consider a two by two density matrix /3 such that is diagonal in the basis f\ £ C 2 , 
fj £ C 2 . Without lost of generality we can consider that 


0 


Pi 0 \ 

0 P2 )' 


PUP2 > 0. 

We will describe initially in coordinates some of the definitions which were used before 
on the paper. 

If 


and 


then 


and 


R®S = 


r=( r " 

*12 

V *21 

*22 


*n *12 \ 


O = 

V *21 *22 ) ' 


/ *1 i 5 i 1 

* 11*12 

* 12*11 

* 12*12 

* 11*21 

* 11*22 

* 12*21 

* 12*22 

* 21*11 

* 21*12 

* 22*11 

* 22*12 

V * 21*21 

* 21*22 

* 22*21 

* 22*22 


\ 

/ 


Tr 2 (*®*) 


*ii(*ii + * 22 ) *i 2 (*n+* 22 ) \ 
* 2 i(*n + * 22 ) * 22 (*n+* 22 ) / 


Given 


( TU 

*12 

*13 

*14 \ 

*21 

*22 

*23 

*24 

*31 

*32 

*33 

*34 

V *41 

*42 

*43 

*44 / 


then, in a consistent way we have 



T 11 + T 22 
*31 + *42 


*13 + *24 \ 
*33 + r 44 J 


The action of an operator U on M2 ® M2 in the basis e\ (g) f\, e2 0 /1, e\ ® /2, e2 ® f2 is 
given by a 4 by 4 matrix U denoted by 


/ K 

11 12 
u u 

r/U 

U 12 

U 11 \ 


If 22 

u n 

f / 21 

U 12 

U\l 

^21 

TJ 12 
U 2\ 

r/U 

U 22 

Ug 

V u 2 \ 

II 22 

U 2\ 

7/21 

U 22 

/ 

(Hi 

TJ 21 

^11 

r/U 

U 21 

M \ 

U 12 
u \\ 

u 22 

U 11 

7/ 12 
U 21 

7/ 22 
U 21 

u\\ 

IJ 21 

U 12 

7/H 

U 22 

^22 

V u\ 2 

TJ 22 

U 12 

TJ 12 
U 22 

t/ll / 
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If U is unitary then UU* = /. This relation implies the following set of equations: 

1) K ufx + U 1 \ 2 u^+ u\\ug+ u\ 2 W = 1 , 

2 ) U\l Uft + U l \ 2 U 22 + U\\uJl+ U\ 2 Ufi = 0, 

3) ull Uf + U l \ 2 uJ 2 + u\\ug + u\lug = o, 

4 ) U\l Wx + K 2 Wx+ U\\W + u\lu% = o, 

5 ) u 2 ' u \1 + u\\ W+ullWi + uflW = ^ 

6) U 2 , 1 Uft + ulfu 22 + ull Wl+U\\uf 2 =\, 

7) Wx + U%Utf + ull Ull + UffUft = 0, 

8 ) uH Wx + UnWx + ull W + ufJJi = 0 , 

9) U\\UU + UUUl 2 + Ull Wi + Uh 2 UJ! = 0, 

10) uffuft + ulluf +ull ull + uim\ = o, 

11 ) u \I W + ul 2 ull + ull ull + ul 2 2 uft= 1 , 

12) u\\W+ uft ufx + ull ull + ull u£ = 0, 

13) ull Wx + ufx Ull + ull W + uftuft = 0. 

14) U 2 \ Wx + ul 2 Uft + ull Wi + uftUft = 0. 

15) u 2 \ Oil + ull Wx + ull Uli + U%U? 2 = 0. 

16) ul\ Wx + ufx ufx + ull W+ufu£= 1 • 

Equation 2) is equivalent to 5), equation 12) equivalent to 13), equation 8) equivalent to 
15), equation 3) equivalent to 9), equation 7) equivalent to 10) and equation 4) equivalent 
to 14). Then we have 6 free parameters for the coefficients of U. 

Using the entries U l r J s we considered above we define 

L(Q)=Pxt 

i=i 

2 

pi E 

i=i 


( ^xx 

Un ' 

\ 0 ( 

' < 

U'h \ 

Wxi 

Uii , 

j Q ( 

, U‘lx 

Uli ) 


( ull 

Ufx 

\ o ( 

' ull 

U?2 \ 

V U ?2 

Uf . 

) Q \ 

, Ufx 

uf ) 
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We can consider an auxiliary L (/ - and express 

m = i(VFi(u n T)Q(V^u il ) + i(Vp-2(u a y)Q(VT2U a ) = 

i'=l i=l 

j^L* l QL n + ^L* 2 QL i2 = £ L*jQLij. 
i=l i=l i.;=l 

From the fact that UU* = I it follows (after a long computation) that 

L{t)=I. 

Note that L preserve the cone of positive matrices. 

Using the entries if‘ r { described above we denote 


£(0) = pi£ 

i— 1 


TJ n TP 1 
U 11 u \2 

TP 1 TP 1 

u 21 u 22 


0 


IT n TJ n 
U 11 U 2l 

TJ n n il 
u 12 u 22 


2 

P2 £ 
1=1 


Tji 2 Tji2 
U \\ U \2 
tji2 t ti2 
u 2\ Ur > 


J 22 


( jri2 t ji2 \ 2 

Q % % ) = Y LjjQL*:. 

I u l2 m 2 1 ^ 11 


'12 u 22 / 1 . 7=1 
One can also show that L(0) = Tr 2 [U (0® /3) U*] (see |3l). 

The first expression is the Kraus decomposition and the second the Stinespring dilation. 

Moreover L preserve density matrices. This is proved in the appendix but we can present 
here another way to get that. If 0 is a density matrix, then 

7>(L(Q)) = Tr{ £ LijQL*j) = £ Tr{L tJ QL*j) = £ MQL*^) = 


1-7=1 


i-7=l 


i-7=l 


Tr( £ QL^Lij) = Tr(Q £ L*^) = 7>(0) = 1 


1-7=1 


1-7=1 


We denote 


Then, 


0 = 


011 012 
021 022 


u ij 0 {u ij y = 


Tj’J TpJ 
u \\ u 12 

if J T/‘j 
u 2l u 22 


011 012 
021 022 


£n U 2 i) = 

u‘A ug ) 


( U^WxQu +U i2 2 Q 2l ) + uUu\[QnfU^Q 2 2) 

V ^ 1 (^ 1011 +^ 021 ) + ^(^012 + t 4022 ) 

We have to compute 


^(£/n0n+t/i202i) +U&(U l AQn + U&Q 2 2) 
uUKQn + U&Q21) + U&iU&Qn + U&Qzi) 


L{Q) = Pl [U u Q(U n Y + U 21 Q(U 21 )*} +p 2 [U 12 Q(U l 2 y+U 22 Q(U 22 y]. 
The coordinate a\ \ of L(Q) is 

Pi [Uti(u!lQn+ullQ2i) + Uti(U 1 1 lQn + u!lQ 2 2)]+ 

Pi [U 2 1 {U\\Qu +U\\Q 2 i) + U\ 2 012 T Q22) ]T 

P2 [Uli(U{\Q\\ +U\ 2 Q 2 i) + U\ 2 (JJ\\Q \2 T Uuf 022) ]t 
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p 2 [UffiuftQn + UftQn) + U\l{U\lQ X 2 + UnQ 22 ))- (2) 


The coordinate a 12 is 

Pi {^(uliQn+u{lQ 2l ) + u£(uilQn + u"Q 22 )]+ 

Pi \U 2 l{UuQn + U\\Q 2 i) + U 22 {U\\Qi2 + U\ 2 Q22)] + 

P2 [U 2 i(uliQn + t/i202i) + U2 2 (ul{Qi2 + u} 2 Q 22 )]+ 

P2 [U%(uftQi 1 + U\ 2 2 Q 2X ) + U%(uf?Q n + U%Q 22 ) ]. (3) 

We will consider a parametrization of the density matrices taking Qn = 1 — Q 22 and 
2i2 = 221- 

The variable £>11 is positive in the real line and smaller than one. Indeed, by positivity 
of £>, we have 0 < 2 11 £>22 = 2n(l - 2n) = 2n ~ Qu_ 

2i2 is in C = R 2 but satisfying 2n(l — 2n) - 2n2i2 > 0 because we are interested 
in density matrices which are positive operators. 

The numbers p\ and p 2 are fixed. Consider the function G such that 

2(2ll>2l2) = 

( P1 [U"(ullQn + U\ l 2 Qn) + U$(u\{Qn + u£(l - 2n))]+ 

Pi [WiiutlQu + u\\ Qu) + u£(uflQn + U$( 1 - 2n))]+ 

P2 [Utf(UftQu + uljQu) + uE(KQi 2 + u! 2 2 (1 - 2n))] + 

P 2 [Uff(uffQn + u\%Qn) + U 22 (U^Qn + u\\{ 1 - 2n))] , 
Pi\^(ullQn+u{lQ^)+u£(ullQi 2 + uti(l-Qn))]+ 

Pi [U^iUuQu + UilQn) + uM{UuQn + Ui 2 (l-Qn)))+ 

P2 [W.iuHQn + UnQn) + U%(U fiQn + u}H 1 - 2n))]+ 

P2 [Ufi(UffQn+U%Q£) +Cf(t/n2i2 + t/f2 2 (l-2n))] ) 

When there is a unique fixed point for G? 


Example: Suppose!/ =e I P <yX ® cr ' =cos(/3) (/<£>/) + i sin(/3) (cr v (g) <y x ). In this case 


/ cos/3 
0 
0 

y i sin p 


0 0 i sin /3 \ 

cos /3 i sin P 0 

i sin P cos p 0 

0 0 cos/3 J 


Therefore, 


G(Qn,Qn) = 

((pi-piQn +P2-P2Q11), 

Pi (cos/3) 2 2i2 +Pi (sin/3) 2 2i2 +P 2 (sin/3) 2 2i2 +P 2 (cos/3) 2 Qi2 ) = 

( l-2n, pi (cos/3) 2 2i2 +Pi (smp) 2 Qn + p 2 (sin/3) 2 ££2 + p 2 ( cosp) 2 Q l2 ) 

One can easily see that given any aglwe have that 2n = 1 / 2, and Q\i = a determine 
a fixed point for G. In order the fixed point matrix to be positive we need that —1/2 < a < 
1 / 2 . 

In this case the fixed point is not unique. 
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II 


It is more convenient to express G in terms of the variables Q\\ £ [0,1], and (a,b) £ 
R 2 , where Q \2 = a + bi. As these parameters describes density matrices there are some 
restrictions: 1/4 > <2n(l — Q\\) > ( a 2 +b 2 ) and 1 >Qn >0 

We denote by Re(z) the real part of the complex number z and by lm(z) its imaginary 
part. 

In this case we get 

G(Qn,a,b) = 

(2n«i + Pi + (an + 012)0 + i (an — a\ 2 )b, 

Re(Qn^ 2 + P 2 + (021 + 0 . 22 ) a + 1(021 — Cl 22 )b), 

Im(Q u a 2 + + (021 + 022)0 + 1(021 -a 22 )b)). 

where 

at = Pi (uUull - Ugutj + Uf[u 2 } - } + 

P2 [uft ull - ufiuli + WiUu - u%uli} , 


Pi = PiWnUn + UuUn ] + P 2 [U\ 2 ull + U$U%] , 

«2 = Pi\uliu I 1 / - W 2 u\\ + u£ull - ] + 

P2[ujfutf - u£u!l + W x u\\~u%u%}, 


p 2 = P 1 [ U22 U X 2 + ] +P2 [ugug + u%ug \, 


r21rr21 


12rrl2 _ 


t22tj22 1 


u=pi[ullu{{ + u%ul{\+ Pl [ullu{l + u%uft\ 


an — pi [ U{{ U H + U xx Uf 2 ]+P 2 [U\{ U12 + U n U ?2 ] 


«2i = pi [u£u{\ + ul\ul\ ] +p 2 [ugutf + uguft} 


r21rr21 


12rrl2 


r22rr22l 


a 22 = Pi [^21 ^12 + ^21^12 ] +P2 [UjlUn + U 2x Uf 2 ] > 
a i is a real number. As 4 > takes density matrices to density matrices we have that / 3 j is 
also real. 

Note that | ctj | < 1 and 1 > / 3 i > 0 . 

It is easy to see from the above equations that (on +012) and i (an — 012) are both real 
numbers. 

We are not able to say the same for (021 + 022) a or i (021 — a 22 ) b . 

In order to find the fixed point we have to solve 

Giiai +J3t + (an+ 012 ) 0 + 1 (an — 012 ) & = <2n 

Qna 2 + p 2 + (021 +a 2 2 )a + i(a 2i -a 22 )b =a + bi, 
which means in matrix form 

f («i ~ 1) an+ai 2 i(a\\—a\ 2 ) \ J ^' | _ / —/3i \ 

\ «2 021+022—1 1(021— 022— 1 ) )\b) V ~P 2 J 

We are interested in real solutions Q\\ ,a,b. 
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In the case of the example mentioned above one can show that cq = 1 and Oo = 0 which 
means that in the expressions above we get a set of two equation in two variables a,b. 
Remember that we are interested in matrices such that 1/4 > Qn(f — Qn) > (a 2 +b 2 ). 
Notice that 0 < Qn < 1. As <J> takes density matrices to density matrices there is a fixed 
point for G by the Brower fixed point theorem. The main question is the conditions on U 
and (3 such that the fixed point is unique. 

If there is a solution ( Qn,a,b ) ^ (0,0,0) in R 3 to the equations 


2n(ai - 1) + («n +a 12 )a+ i(an -a n )b = 0 
Qndi + (ci2i+C 122 —3)a+i(ci2i — ci22~ 3)b =0, (4) 

then, the fixed point is not unique. The condition is necessary and sufficient. 

A necessary condition for the fixed point to be unique is to be nonull the determinant of 
the operator 

K= ( +an ; (an — an) \ 

\ «21 4“ «22 1 fl22~ 1) ) 

Notice that if (zi,Z 2 ) satisfies K{z\ 7 Z 2 ) = (0,0), then ^ is real (because an +«i 2 and 
/(an — a n) are real). From this follows that there exist a solution ( a,b) £ R 2 in the kernel 
of K. In this case (0 ,a,b) is a nontrivial solution of 0. 

The condition det K ^ 0 is an open and dense property on the unitary matrices U. 
Indeed, there are six free parameters on the coefficients U' r ] s . Consider an initial unitary 
operator U. One can fix 5 of them and move a little bit the last one. This will change U 
and will move the determinant of Kjj in such way that can avoid the value 0 for some small 
perturbation of the initial U. 

Suppose U satisfies such property Det U 0. For each real value Q\\ we get a different 
(ag n ,/7gu) which is a solution of K{a 1 b) = (— Qu (ai — 1), — Qn 0 C 2 )- 

In this way we get an infinite number of solutions {Qn,aQ n ,bQ n ) S K x C 2 to 0. 

a 2 is not real. 

But, we need solutions on M 3 . Denote by S = Sy the linear subspace of vectors in C 2 
of the form p(cq — 1, OC 2 ), where p is complex. 


Lemma 3. For an open and dense set of unitary JJ we get that K 1 (S) fi R 2 = {(0,0)}. 
For such U, suppose {Qn,a,b) satisfies equation (0, then the non-trivial solutions ( a,b) 

°f 

K(a 7 b) = (-Qn(ai - 1),-Qna 2 ) 

are not in R 2 . 


Proof: Suppose 1 f //' = a + J3 / = z° = Zg. Note that for a generic U we have that 

d2 -fi 0. 

We denote C\\ — a\ 1 -\-a 12 , G\ 2 — / (flu — fli 2 )> On = ^ 2 i~ffl 22 — 1 and finally C 22 = 
/(a 21 —022 — 1)- 

Suppose ( Qn,a,b) € M 3 satisfies equation 0. We know that generically on U the value 
Qn is not zero. 

For each C (/ - we denote C,j = Cjj + Cfj /, where i,j= 1,2. 

If K(d,b) = (-Qn(ai - l),-<2n a 2 ), then 

Cnd + C2i^ = z° (C2\d + C22b) = (a + / 3 /) (C2i« + C22^)- 

In this case 

Cnd + C2i^ = (aC^a — PC^a — fiC^b — dC22b)+ 
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i (/3C 21 a + ac\ x a + aC\ 2 b — fiC 22 b). 

If a and b are real, then, as C\ \ and C 22 are real, then 

(/3 C l 2l + aC 21 ) a + ( aC l 22 - j3 C\ 2 )b = 0. (5) 

Moreover, 

(aC l 2l - J3C 2 2 ! - C n )a - (j3Ci 2 - aC 2 2 2 - C 21 ) £ = 0 (6) 

If 


Det 


fiC 21 H - 

«C 21 -j3C 21 — Ci 1 


aC, 1 , - 6C 2 2 
j3Ci-aC 2 -C 21 


± 0, 


then just the trivial solution (0,0) satisfies ([5]) and ©. 

The above determinant is non zero in an open and dense set of (/. 
Then, the solution (Q\ 1 , a,b) £ ffi 3 of ((4]) have to be trivial. 


□ 


Under this two assumptions on U (which are open and dense) the fixed point for G is 
unique. Then, it follows that the density matrix Q = (9<t> which is invariant for <t> is unique. 
Given an initial Qq any convergent subsequence O' 1 '- ((jo), K —> °° will converge to the fixed 
point (because is unique). 

As 

G(Qu,a,b) = 

(filial + Pi + (#11 +^12)0 + *(<211 — ^12) b, 


Re ( Qi 1 a 2 + Pi + {a 2 i + 022 ) a+i(a 2 \-a 2 2 )b), 
Im(Q n a 2 + 1 3 2 + (<221 + a 22 )a + i(a 2 i -a 22 )b)), 

one can find the explicit solution 


2<t> = 


Q 11 a + bi 

a — bi 1 — Qn 


by solving the linear problem G(Qn,a,b ) = (Qn,a,b). 


4. Appendix 

Lemma 4. Given A, B £ Jf(V), then Tr(A®B) = Tr(Tr 2 (A®B )). Moreover, Tr{ Tn(T )) = 
Tr(T), for all T £ «£?(V (g> V). 

Proof: 

Indeed, 


7>(A ®B) = Tr(A)Tr(B) = Tr( Tr(A)B) = Tr( Tr 2 (A®B)). 


□ 


Lemma 5. Given T £ Jz?(V ® V), 

a) if T is selfadjoint, then, Tr 2 is also selfadjoint, 

b) moreover, ifT is also positive semidefinite then Tr 2 {T) is semidefinite. 
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Proof: 

a) If T is selfadjoint, then, = t k nj. This implies that T.j f ijkj = L/ hjij ■ Therefore, 
7>2 is selfadjoint. 

b) If T is postive semidefinite, then < T(x®x') , x ® x J >> 0, for all x,x! £ V. In 
particular, < T(x® e q ), x® e q >> 0, for all x = ci e\ + .... +c„e n £V and 1 <q<n. 

As T(x®e q ) ='Ltij k iL ik {x)®L j i{e q ) = 'Lti jkq c k (e i ®ej )., then 

< T(x®e q ),x®e q > = Y, t iqkqC k cl, q= 1,2, 
i,k 

From this follows that Li,k,qt iq kqCkCi = Li.k CLqtiqkq) C k c]>0. 

Then, < Tr 2 (T) ( x),x >> 0. 

□ 

Note that the analogous property for positive definite T is also true. 

Lemma 6. If A £ T, then Oy (A ) £ T, for all U £ f/. 

Proof: 

As A and /3 are selfadjoint and positive semidefinite the same is true for ,4 0/3, Then, 
the same is true for f/(A® fi)U*. From Lemma|5]we get that <T>{/(A) = Tr 2 (U{A®fi)U*) 
is selfadjoint. 

By Lemma[4]7>(<I>j/(A)) = Tr(U(A®p)U*) = Tr(A®b) = Tr{A)TR(B) = 1. 

□ 
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